The ordering temperature and critical exponents of the binomial Ising spin glass in 

dimension three 
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We compare numerical estimates from different sources for the ordering temperature Tg and 
the critical exponents of the Ising spin glass in dimension three with binomial (±J) interactions. 
Corrections to finite size scaling turn out to be important especially for parameters such as the 
Binder cumulant. For non-equilibrium parameters it is easier to approach the large size limit and 
to allow for corrections to scaling. Relying principally on such data, a crossing point defines the 
freezing temperature Tg; the possibility that the ordering temperature is zero can definitively be 
excluded. We estimate an ordering temperature Tg = I.f95(f5), with associated estimates of the 
critical exponents for which corrections to finite size scaling are well under control. Among the 
parameters evaluated is the leading dynamic correction to scaling exponent w. 



INTRODUCTION 

The Edwards- Anderson Ising spin glass (ISG) has 
been studied intensively for over twenty five years. In 
this model Ising spins on a regular lattice in dimension 
three are coupled by random near neighbour interactions 
101 . Despite its deceptively simple definition, the ba- 
sic physics of this canonical complex system still eludes 
consensus. Recently interest has focused on the correct 
description of the ground state and the low energy ex- 
cited states 0, but finite temperature aspects are also 
highly important, in particular the critical behaviour at 
the freezing temperature Tg. Over the years a number of 
numerical efforts have been consecrated to accurate mea- 
surements near Tg, with the aim first of demonstrating 
that Tg is non-zero, and then of determining the critical 
exponents precisely ^, ^, ^, 0, ^ . The long relaxation 
times and the need to average over large numbers of inde- 
pendent samples renders the task numerically laborious. 
In addition, it has not always been appreciated that cor- 
rections to finite size scaling can be singularly vicious in 
this system. 

Here we will concentrate on the problem of the ac- 
curate and reliable determination of Tg and the critical 
exponents for the ISG with binomial (±J) near neigh- 
bour interactions in a dimension three simple cubic lat- 
tice {L X L X L), 
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with Ji j chosen at random as ±1. As this system is 
the dropsophilia of spin glass studies, it is important to 
establish as precisely and reliably as possible the basic 
parameters of its critical behaviour. We will show that a 



parameter set can be obtained consistent with most pub- 
lished data from different complementary approaches, if 
corrections to finite size scaling are properly taken into 
account. 

We will principally refer to the simulations of refer- 
ences [@, H, 0, H ; the estimates drawn from the analy- 
ses of the numerical data are shown in Table I. We will 
comment on these data as we go along. The main nu- 
merical contribution of the present work is to present 
non-equilibrium data of higher statistical accuracy than 
before, and to allow for corrections to finite size scaling 
in the analysis of these results also. 

As a general rule, the apparent exponent values are 
strongly correlated with the value assumed for the or- 
dering temperature, so once Tg is well determined it is 
relatively easy to obtain good values for the exponents. 
On the contrary, if the wrong Tg is used, the estimates 
of the exponents will be biased. Finite size scaling meth- 
ods have been widely used, but corrections to finite size 
scaling have often been overlooked, which is dangerous. 
In Table I the estimates of the ordering temperatures 
and of the exponents from different publications should 
not be considered as independent of each other, but each 
parameter set should be taken as a whole. 

The autocorrelation function at time t is defined by 



qit) SmMt) > 



(2) 



where the average (< .. >)is over all spins. In practice 
a large number of further averages (denoted by [..])are 
taken over many samples all of size spins with micro- 
scopically different arrangements of the bonds. We will 
use q to mean q(t) at times long enough for the system 
to have entered the equilibrium fluctuation regime. Pe- 
riodic (toroidal) boundary conditions are generally but 
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not always used (see j^). The equilibrium spin glass sus- 
ceptibility is given by 

Xl{T) = L^[< >] (3) 

where the average is taken over time once the sample 
has been annealed to equilibrium either by carrying out 
enough annealing steps or by using the more sophisti- 
cated and efficient parallel tempering method ||^, ^ . The 
non-equilibrium susceptibility xi^) is the spin glass sus- 
ceptibility at time t 

X{t)^L^[<q\t)>] (4) 

after a simple anneal which starts from a totally uncorre- 
lated (infinite temperature) state and lasts a time t . 

BINDER PARAMETER 

The first large scale simulations on this problem were 
those of Ogiclski Q] who did massive dynamic and static 
measurements on samples as large as L = 64. His val- 
ues were considered authoritative until finite size scaling 
techniques were later intensively exploited ||, ||, ^, |j. 
The standard scaling method to estimate Tg from finite 
size data in ISGs is through the use of the dimensionless 
Binder parameter, 

5l(T) = [?,- <q^> / <q^ >2]/2 (5) 

where < q^ > and < q"^ > are moments of the fluctu- 
ations of the autocorrelation parameter q{t) in thermal 
equilibrium at temperature T. The averages are taken 
first over time and then over samples. In absence of cor- 
rections to finite size scaling, the family of curves gL{T) 
for different sizes L must all intersect precisely at for 
T = Tg. For the 3d binomial ISG, the gL{T) curves come 
together near T = 1.20 but very good statistics are 
needed to see intersections |5[ ^. Precisely because the 
gL{T) curves lie very close together, even small correc- 
tions to finite size scaling can have a drastic effect on the 
crossing point temperatures. As corrections to finite size 
scaling are already visible for the spin glass susceptibil- 
ity j^, 1^ , corrections can be expected a fortiori for the 
Binder parameter. If uj is the leading correction to scaling 
exponent, at Tg the successive gL values for increasing L 
can be expected to behave as 

gL ^gn[l + const ^L-"^] (6) 

It has been pointed out j|] that in these circumstances 
it is helpful to use the " quotient method" . An elementary 
version of this method is to use " crossing points" . Succes- 
sive intersection temperatures T*(L), corresponding to 
intersections of curves gL{T) and g2L{T), are plotted. If 
the leading correction exponent lo dominates, then T*{L) 
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FIG. 1: A selection of results for the Binder parameter gL{T) 
as function of the temperature T. Data for L = 5, 10, 20 from 
data for L = 3, 4 present work. Error bars are about the 
size of the symbols. 

values should extrapolate to the true infinite L ordering 
temperature Tg as 

T* (L) =Tg + const * £-("+1/-^) (7) 

Extrapolation to infinite L provides a precise esti- 
mate for Tg, if the corrections to finite size scaling are 
fully under control. In the case of the 3d binomial 
ISG, accurate gL{T) data for sizes L = 5, 10,20 for 
L = 6, 8, 12, 16, 24 || and for L = 3, 4 § can be used to 
obtain intersections. It should be noted that for technical 
reasons the boundary conditions used by Ballesteros et al 
Q on a specially built dedicated computer are helicoidal 
while for the other calculations the boundary conditions 
are periodic. A direct comparison of the raw data from 
§ and 1^ indicates that systematic differences in the 
results because of the different boundary conditions are 
negligible for the spin glass susceptibility at L = 20. The 
sequence of Binder parameter curves for L = 3, 4, 5, 6, 8 
1^, 0, 1^ behaves smoothly, although the L = 5 curve is 
the only one in the set to have helicoidal boundary con- 
ditions. A subset of the data from the different sources 
is shown for illustration in figure 1. 

It can be seen that the behaviour of the 5l(T) curves 
and hence of the crossing temperatures varies in a non- 
monotonic manner with L. (We suggest below that this 
may be because uj is high and analytic corrections may be 
non-negligible) . The ui and v values of Ballesteros et al |^ 
correspond to an effective exponent {oj + l/i/) — 1.3(4). 
(Estimates of to vary considerably from author to author 
and we will return to this point later). In figure 2 we 
have plotted T*{L) in the form of equation (2) using the 
central value of the exponent. The pairs of sizes [2L, L] 
we have used are [20,10] and [10,5] where the bound- 
ary conditions are helicoidal, [16, 8], [12, 6], [8, 4] and [6,3] 
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FIG. 2; Quotient plot of crossing points for the Binder pa- 
rameter. Tcross{L) is the temperature at which the curves 
QhiT) and g2L{T) cross. The x axis is L~^ '^ (see text). The 
curve is just a guide to the eye. Tcross{L) should extrapolate 
linearly to Tg for infinite size L as 

where the boundary conditions are periodic. There is 
no very obvious systematic difference between the two 
sets of points with different boundary conditions. Even 
for the larger sizes the crossing temperature is still in- 
creasing fairly rapidly with increasing L ; the very high 
quality data of show a L = 10, 20 crossing already at 
T* = 1.118(10), so the central estimate given by [||, 
Tg ~ 1.11, is definitely too low. The curve in figure 2 
(a two paramereter spline fit) is just a schematic indica- 
tion of one possible extrapolation to to infinite L which 
would appear to lead to a significantly higher estimate for 
Tg than that of g, Tg = 1.138(10). Alternative curves 
could be drawn through the data points, but the Binder 
parameter data (including the high quality results of j|] ) 
point to a higher Tg than their central estimate. The non- 
monotonic behaviour of the crossing points means that 
not only the leading correction term is operative. Thus 
in the case of this particular ISG system the Binder pa- 
rameter method can only give a rough estimate of Tg, 
because the correction terms are badly under control. 

It should be noted that because the behaviour of the 
crossing temperature T*{L) is not monotonic, if only 
measurements of gL{T) for i = 4, 8, and 16 were avail- 
able one would see an excellent unique intersection of the 
three gL{T) curves, but this intersection point is NOT at 
the true infinite size limit temperature Tg. It is an arte- 
fact of the non- monotonic behaviour of T* (L) . 

Binder parameter data can also be directly used to 
estimate the exponent v, as 

%i)Ml/T) cc L^'^ (8) 

at Tg. Using a single gL{T) data set, the value of v esti- 
mated in this way depends strongly on the Tg estimate. 
The variation of the gL{T) slopes with L in the gL data 



of Ballesteros et al provide an apparent v value which is 
about 2.15 if Tg = 1.138 (their preferred value); if Tg is 
higher then the associated estimate for v from the same 
data set is lower. Once again estimates can be affected 
also by correction terms. (Similarly, the apparent 77 from 
spin glass susceptibility data changes rapidly with the 
assumed Tg, Figure 7). 



CORRELATION LENGTH 

Ballesteros et al Q do not in fact base their Tg esti- 
mate on their Binder parameter data, but use an inde- 
pendent criterion : the intersection of curves for the ratio 
of the size dependent correlation length to L. However 
they show data for three sizes only {L = 5, 10, 20) and 
do not discuss possible corrections to finite size scaling 
for this parameter, even for L as low as L = 5. Even 
though their intersection point for these three sizes ap- 
pears to be unique and very clear, the Tg estimate from 
this criterion {Tg = 1.138(10)) is lower than the prefered 
Tg value obtained from other measurements which we 
will discuss below, Tg ~ 1.195(15). Although both series 
of measurements show clear evidence for a well estab- 
lished transition, and though the difference between the 
estimated transition temperatures is only about 4%, it 
would be more satisfactory if complete agreement could 
be reached. 



SERIES 

It is of interest to note that quite independent infor- 
mation from series calculations can throw further light 
on the values of the ordering temperature and the ex- 
ponents. Long series results on the 3d binomial system 
|]lO| provide sets of parameters \Tg,^]a corresponding to 
different approximants. Figure 3. The parameter pair 
[rg,7] representing the true phase transition can be ex- 
pected to lie somewhere within the cloud of approximant 
points. In Figure 3 we represent the approximant points 
from |l^ together with the equivalent points from the dif- 
ferent numerical simulation estimates from Table I (using 
7 = j/(2 — 77)). It can be seen that while the approximant 
cloud has a relatively wide spread in temperature and so 
is not very discriminatory for Tg, the cloud has a fairly 
narrow distribution in 7. The 7 values corresponding to 
the parameter sets of references [||, ^and |^ in Table I 
lie systematically above the series data, while the [T'g,7] 
pairs of and ||^ fall neatly in the centre of the approxi- 
mant cloud. In dimension 4 the agreement between series 
pO| and simulations jlll is excellent, so there seems no 
a priori reason to expect the series results to be strongly 
biased in dimension 3. 
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FIG. 3: The estimated critical exponent 7 plotted against 
the estimated ordering temperature Tg. Squares are values 
obtained for different approximants from series calculations 
[ p^ ; the couple [79,7] corresponding to the true ordering 
temperature and exponent for the system can be expected 
to lie somewhere within the cloud of points indicated by the 
ellipse. Circles are values reported from different simulations; 
from left to right : |,|], §,§,§. 

NON-EQUILIBRIUM METHODS 

We have proposed a novel technique including non- 
equilibrium data for identifying Tg numerically through 
the combination of three simulations |0, |l2| . We will use 
it again in the present context in a slightly modified form. 

First, the non-equilibrium spin glass susceptibility for 
samples quenched to Tg increases with time as t'' where 

h={2- ri)lz (9) 

T] and z being the standard equilibrium static and dy- 
namic exponents |p^ . 

Secondly the autocorrelation function decay q(t) at Tg 
is of the form , if a sample is initially annealed at 
Tg during a long waiting time tw^ under the condition 
that the measuring time for the decay t is much shorter 
than tw ^, 0. The exponent x at Tg is related to the 
equilibrium exponents through x — [d — 2 + Tf) /2z 

Thirdly, the equilibrium spin glass susceptibility at Tg 
increases with sample size as 

Xl/L^ cx L-^ (10) 

The effective exponent h{T) can be measured directly 
as a function of trial temperatures T. It can also be eval- 
uated indirectly and independently through combining 
the effective exponent "qiT) from the finite size scaling 
of the equilibrium susceptibility, with the effective decay 
exponent x{T). We will refer to this derived value as 
h*{T) : 

h*{T) = 2x{T){2 ~ r,{T))/(d - 2 + ry(T)) (11) 



Measurements of each of the three effective parame- 
ters {h{T),x{T),ri(T)) can be made at a series of trial 
temperatures T; consistency dictates that at Tg we must 
have h^h* so when h{T) = h*(T), then T ^Tg. 

Once again, care must be taken to avoid errors arising 
from corrections to finite size scaling. For the measure- 
ments of h{T), large samples can readily be used, so it is 
relatively easy to insure that for the time range used the 
system is far from the saturation limit for that particular 
size, see At short times there can however be finite 
time corrections, the analogue of finite size corrections, as 
the clusters of correlated spins that are gradually build- 
ing up with annealing time t have a t dependent finite 
size, even if the sample size L can treated as effectively 
infinite. Including the finite time correction, we should 
write |l|l 

x{t)=At''[l~Bt-'"/'] (12) 

with h{T) as above, and w the leading dynamic correc- 
tion to scaling exponent (which does not have to be equal 
to the static exponent cj[0). The factors A and B may 
vary with T. An example of data is shown for two alter- 
native plots in figures 4 and 5. (Rather than x{t) we plot 
x{t) — 1 because at infinite T, x{t) = 1- Measurements 
were averaged over runs on 6766 samples of size L = 28 
at temperature T = 1.2. For this size saturation to the 
finite size limited equilibrium susceptibility would only 
occur at i ~ 10*). For this particular temperature the 
best fit gives for the important exponent h{T) a value of 
0.390(5). Our results for h{T) at different temperatures, 
figure 7, are in good agreement with those of |l^, |l8| , but 
have higher statistical accuracy and allow for the correc- 
tion to scaling which was not considered in the earlier 
work. We have checked that our estimates for h{T) do 
not change when we use sample sizes L from 20 to 28; 
this means that we should be in the limit such that the 
h{T) values represent the infinite size behaviour. 

Estimates of w/z from the analyses of the non- 
equilibrium susceptibility data have statistical fluctua- 
tions but do not show any systematic trend with T over 
the fairly narrow range of T values near Tg that we have 
explored, so our estimate of w/z is not sensitive to the 
exact choice for Tg. The mean value is w/z = 0.48(10). 
As a number of independent measurements have shown 
that the dynamic exponent z is near 6 ( 0, |l^ and our 
estimate given below at the end of the analysis), we can 
deduce from w/z that the dynamic correction to scaling 
exponent w ~ 2.9(6). This turns out to be very close 
to the static exponent value lo ^ 2.8 estimated from an 
analysis of the size dependence of the spin glass suscep- 
tibility x{L) 0. It is distinctly higher than the estimate 
uj = 0.8(4) given by (An attempt to fit the q'^{t) data 
with a correction exponent w/z of the order of 1/6 fails 
completely) . 

There is an independent check on the expected size of 
UJ. Estimates of uj can be obtained from series expansion 



5 




10 100 1000 10000 

t 



FIG. 4: An example of the extraction of the Huse expo- 
nent h{T) @ from non-equilibrium susceptibility x(i) data. 
(x(i) — 1) is plotted against time t on a log/log scale. The 
data in this particular plot correspond to T = 1.2 and are an 
average over 6766 independent samples. Data : squares, fit : 
circles 
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FIG. 5: The same data as in Figure 4, but presented in the 
form of a (x(i) — 1)/^'' against log{t) plot with h{T) chosen as 
the best overall fit value. The values of h,w/z,A and B are 
adjusted to give the best fit to (xW-l)A'' = A(H-Bt""'/^). 
Data : squares, fit : line 

data through uj = Ai/u ||l9| where Ai the appropriate 
correction to scaling parameter for the series work. In 
dimension 4 a reliable value ^ 3.0 was obtained [ p^ ; 
in dimension 3 the estimate by the same authors is ~ 3.0 
with rather lower accuracy. In the sequence of Ising ferro- 
magnets and for the sequence of percolation models 
pl[ , Lo increases regularly with d„ — d, where d„ is the 
upper critical dimension (4 for Ising ferromagnets, 6 for 
the percolation model as for ISGs). The high value of oj, 
near 3 in dimension 3 from the ISG series work, is consis- 
tent with the value from the present simulation analysis 
and is also consistent with there also being a general up- 
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FIG. 6: The autocorrelation function initial decay exponent 
x{T)- data from 

ward trend for uj with increasing d„ — d in ISGs. A much 
lower estimate for lo Q seems inconsistent with the series 
work. 

A large correction exponent means that the correction 
factors tend towards 1 quite rapidly as L or t increase, 
so the large size or long time limiting behaviour can be 
reached fairly readily. However for large to it may not be 
safe to ignore analytic corrections to scaling which 
could perhaps explain the non-monotonic behaviour of 
T£ from the Binder parameter data seen above, figure 2. 

For the estimates of x(T), the measurement runs for 
q{t) must be preceded by long annealing runs. It turns 
out that after an appropriate anneal the algebraic form 
with constant x{T) sets in from rather short times (a few 
Monte Carlo steps) , and the measurements are self aver- 
aging jj, |l4j . Deviations from the asymptotic behaviour 
are so small that no meaningful analysis in terms of cor- 
rections to scaling could be carried out (see for instance 
the q{t) data shown in jj]). As a consequence it is rela- 
tively easy to obtain robust and accurate values of x(T) 
on large samples. The data of references jj] and are 
accurate and in excellent agreement with each other, fig- 
ure 6. Our own data are in full agreement with these val- 
ues so we have relied on the published measurements. In 
ISGs the relaxation function continues to be algebraic for 
temperatures below Tg, so there is no problem in defin- 
ing x{T), while above Tg an additional relaxation term 
with a cutoff function form appears. Ogielski introduced 
phenomenologically the stretched exponential decay in 
this context to parametrise the cutoff function, and it 
has been frequently used since. For present purposes, at 
temperatures where the pure algebraic decay no longer 
holds, one is clearly above Tg. 

Again, as for h(T), x{T) can be measured on almost 
arbitrarily large samples. This is an important advantage 
of methods where achieving strict thermal equilibrium is 
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FIG. 7: The effective value of the equihbrium exponent r]{T), 
eqn (3), from finite size scahng data 



not obligatory. 

Numerically the most demanding of the three measure- 
ments is that of the size dependent equilibrium spin glass 
susceptibility x(L,T). We have relied on the excellent 
data of H], which extend up to sizes L from 16 to 24 
depending on the temperature. The method of allowing 
for corrections to finite size scaling is outlined in Q . Es- 
sentially, a reliable estimate of the effective temperature 
dependent exponent rjiT) can be obtained from a fit to 
x{L,T) including the correction factor 



(13) 



with an effective exponent i^lT) at and below Tg and a 
temperature independent uj. The correction term turns 
out to be more important than the statistical error only 
for L < 8. (We have checked that the raw equilibrium 
xiL,T) data exhibited by ||^ are in excellent agreement 
with those of [^). The x(L,T) results can be taken to 
provide good infinite size limit estimates for r]{T), fig- 
ure 7, with individual errors of about ±0.03 including 
statistical and systematic contributions. Above Tg the 
x{L,T)/L'^ curves start to bend down with large L, as 
would be expected from the general scaling form. 

The results for h{T) and h*{T) as defined above are 
shown in figure 8. There is a clean intersection at 
T — 1.195(15) which we can take to be a reliable esti- 
mate for Tg, essentially free of corrections to finite size 
scaling because in each data set care has been taken to 
use as large samples as possible, and to allow for cor- 
rections to finite size scaling whenever appropriate. The 
clean intersection should set at rest any lingering suspi- 
cion that there is not a true transition in this spin 
glass at a critical temperature Tg close to that estimated 
by Ogielski [Q many years ago. 

We obtain the corresponding exponents from an anal- 
ysis of the various data, with Tg in hand, Table I. Con- 
sistent and precise values are obtained using finite size 



FIG. 8: The Huse exponent h{T) measured directly as in Fig- 
ures 4 and 5 (squares), and the same exponent h*{T) mea- 
sured via the parameters x{T),ri{T), eqn (4) (triangles). For 
consistency, the intersection of the two curves must occur at 
the ordering temperature Tg. 



TABLE I: Estimates of Tg, v, r) and 7 for the 3d ISG with 
binomial interactions using different techniques (see text) 



Reference 



Tg 



Ogielski [| 
Singh- Chakravarty |1 
Kawashima- Young 
Palassini-Caracciolo 
Mari-Campbell 
Ballesteros et al 
Present analysis 



1.175(25) 
1.18(8) 
1.11(4) 

1.156(15) 
1.19(1) 



1.3(1) -0.22(5) 



1.7(3) 
1.8(2) 
1.33(5) 



1.138(10) 2.15(15) 
1.195(15) 1.35(10) 



-0.35(5) 
-0.26(4) 
-0.22(2) 
-0.337(15) 



2.9(3) 
2.8(4) 
4.0(8) 
4.1(5) 
2.95(15) 
5.0(4) 



-0.225(25) 2.95(30) 



scaling for x{L, T) [|j with correction terms [Q, and from 
the large sample data shown by Ogielski Q]. The value 
estimated for z directly from the critical values of h and 
X is 5.65(15) which can be compared with 6.0(8) from 
Q. The major source of error for all the exponents (as 
reflected in the error bars) is the estimated uncertainty 
in Tg. 

Finally, if we plot 7/(d) as a function of dimension d 
we obtain the curve shown in Figure 9, where we have 
included the present result (Table 1) for dimension 3, 
and values from 1 23 1 , ^ , |jl^, ^ and the mean field 
value(77 = 0) for dimensions 2,4,5 and 6 respectively. 
(We can note that rj remains well defined even in dimen- 
sion 2 where Tg tends to zero ||2^ ) . The curve is a guide 
for the eye. The overall shape of ri{d) is strikingly similar 
to that for the analagous series of r]{d) values in the per- 
colation model |2^. The present estimate for ri{d = 3) in- 
trapolates smoothly into the sequence of values obtained 
independently for the other dimensions. 
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FIG. 9: The exponent 77(d) as a function of dimension d with 
values taken from references as indicated in the text. The 
curve is a guide for the eye. 



CONCLUSION 

In conclusion, we have once again considered the prob- 
lem of obtaining reliable estimates for the ordering tem- 
perature of the binomial interaction ISG in dimension 
three and for the associated critical exponents. It is very 
important to allow for corrections to finite size scaling 
for this particular system. It is clear from figures 1 and 
2 that Binder parameter crossing points in this ISG are 
strongly affected by substantial corrections to finite size 
scaling even up to and beyond sizes as large as L = 20. 
For the spin glass susceptibility x, corrections to scaling 
are also present but their influence virtually dies out by 
L 8. From non-equilibrium susceptibility data we ob- 
tain an estimate for the value of the dynamic correction 
to scaling exponent which is large, w = 2.9(6), and very 
similar to a previous estimate of the static correction to 
scaling exponent uj Q. 

From an analysis of non-equilibrium data on large sam- 
ples together with equilibrium susceptibilites |@, |l^ , our 
estimates for Tg = 1.195(15), and the associated static 
critical exponents, are in excellent agreement with those 
of Ogielski who made dynamic and static measure- 
ments on even larger samples (up io L — 64). They are 
close to those of Palassini and Caracciolo , who care- 
fully extrapolated to infinite size, explicitly allowing for 
corrections to finite size scaling. They are consistent with 
high quality series results [ p^ , which come from a method 
complementary to and independent of the Monte Carlo 
numerical simulations. When Binder parameter data are 
extrapolated to large L, the Tg estimate obtained is im- 
precise because of uncertainties in the corrections to finite 
size scaling, but within the errors it is consistent with the 
estimates using the other techniques. There is however a 



puzzling residual disagreement with the work of Balles- 
teros et al |^ concerning the precise value of the critical 
temperature Tg. 

We gratefully acknowledge the use of the computing 
facihtics of IDRIS. 
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